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Abstract 

PT— symmetric potentials V(x) — — x 4 + iB x 3 + C x 2 + \Dx + iF/x + G/x 2 are quasi-exactly 
solvable, i.e., a specific choice of a small G = G^ ES ^ — integer/ 4 is known to lead to wave functions 
^ ES \x) in closed form at certain charges F = F^ ES ^ and energies E = E iQES K The existence 
of an alternative, simpler and non- numerical version of such a construction is announced here in 
the new dynamical regime of very large G^ QES "> — v 00. 
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1 Introduction 



Eight years ago Bender and Boettcher [1] considered a specific VT— symmetric "asymptotically 
repulsive" oscillator 

H = p 2 + V (BB \x), V (BB \x) = -x 4 + 2iax 3 + (a 2 - 2b) x 2 + 2i{ab - N) x (1) 

and conjectured and numerically verified that it possesses the real and discrete spectrum in cer- 
tain intervals of couplings a and b (note that while V denotes the operator of parity, the complex 
conjugation T mimics time reversal so that a and b must be chosen real). These authors empha- 
sized that their non-Hermitian model may be considered, in a way, a "nearest neighbor" of the 
harmonic oscillator as it exhibits, in a sharp contrast to its undeservedly more popular Hcrmitian 
and asymptotically growing +x 4 alternative [2], the exceptional quasi-exact solvability (QES, [3]). 

By definition the latter feature means that in a way which parallels harmonic-oscillator wave 
functions ip(x) ~ exp(— x 2 /2) x a polynomial, a part of the set of the bound states generated by 
the Hamiltonian (1) remains elementary, 

N 

^ BB \ X ) = e -i* 3 /3-az 2 /2-i6* ^ ^ ( 2 ) 

This observation acquires a particular appeal in the light of the recent increase of interest in 
the possible applications of non-Hcrmitian models in quantum optics [4] and in the analyses of 
quantum chaos [5] as well as in various innovations of supersymmetric [6] , magnetohydrodynamical 
[7] or particle-physics [8] models. During the recent quick development of the related theory of 
VT— symmetric models [9] - [11] it has been, moreover, revealed that their quantum bound states 
may be assigned the standard probabilistic interpretation, provided only that one re-defines the 
scalar product in Hilbcrt space, {-\-)uirac — > (t) adapted- For this purpose one only has to introduce 
an unusual, Hamiltonian-dependent metric operator 6 ^ I in a way which proved productive in 
nuclear physics [12], 

(-01^2) Dvrac > {ll>M 2) adapted = ("01 \Q\lp2) Dirac ■ (3) 

It is now agreed [13] that the VT— symmetric Hamiltonians - symmetric) ma y ^ e use( j as 

phcnomenological models whenever we succeed in an explicit construction of the metric operator 
9 = e(H) 1 . 

The latter observations enhance the importance of the Bender's and Boettcher's partially solv- 
able two-parametric model (1) as well as of its straightforward three-parametric "charged" and 
"spiked" generalization, with a Coulomb and centrifugal force added in ref. [14], 

V(x) = -x 4 + i B x 3 + C x 2 + i D x + i F/x + G/x 2 . (4) 

Unfortunately, the phenomcnological applicability of both these models proved unexpectedly hin- 
dered by the computational difficulties arising during the explicit construction of their exact bound 
states (cf. section 2 for a brief review) . In a reaction to such a contradictory situation we returned 

1 See, e.g., the September special issue of Czech. J. Phys. 55 (2005), pp. 1045 - 1192 for more details. 
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to this class of models once more. We revealed and report here a significant simplification of the 
QES construction which emerges in a strongly spiked limit, i.e., for very large couplings G — > oo. 

In full detail our observations will be described in section 3 and summarized in section 4 em- 
phasizing that the new dynamical regime is complementary to the two-parametric option (1) of 
ref. [1] with vanishing G. Our new construction may even be considered simpler since it leads to 
the determination of the negative-quartic VT— symmetric QES bound states in terms of closed and 
compact formulae. 

2 Quartic models and their quasi-exact solvability 

2.1 A broader family of the next-to-harmonic models? 

Among all the available exactly solvable versions of Schrodinger equation H \tp) = E in Quantum 
Mechanics, an undoubtedly exceptional position belongs to the harmonic oscillator, the Hamiltonian 
of which preserves the same differential-operator form in both the x— and p— representations [15]. 
Although such a curious "Fourier-transformation-symmetry" property of H^ H °^ = p 2 + x 2 does 
not survive the transition to the "next", quartic anharmonic oscillators, it still may play a role in 
their pcrturbative [16] or continued-fraction [17] description. Moreover, an unexpected role of the 
Fourier-transformation partnership between two different quartic oscillators has been revealed by 
Buslaev and Grecchi who succeeded in proving a strict isospectrality between certain two "next-to- 
harmonic" quartic-oscillator models H^ Hermitian 1 and H^ vr -^ mmetri ^ (cf. [18]). 

The subsequent increase of interest in VT— symmetry in Quantum Mechanics [9] climaxed re- 
cently, in the specific quartic-oscillator context, with the paper [19] where, for a sample choice of 
the negative-quartic H^ VT ~ s v mmetrlc ) ^ — x 4 , an explicit construction of the metric 9 has been 
presented as performed without ad hoc tricks and starting simply from the first principles. The re- 
lated Buslaev's and Grecchi's results are recollected there so that, in some sense, the circle is closed 
and the picture seems completed. Yet, the description of another, viz., QES harmonic-oscillator-likc 
property of models H^ VT ~ symmetrlc "> deserves an independent completion. 

In a way indicated by Buslaev and Grecchi ([18], cf. also [20]) and re-emphasized, e.g., by Dorey 
et al [21], our understanding of the various aspects of VT— symmetry may be made simpler rather 
than more complicated by an introduction of the angular momentum L in our ordinary Schrodinger 
equation, 

tP(x) = Eip{x) . (5) 
Traditionally one abbreviates 

L=l(d-3), l + l(d-3), 2 + ±(d-3), ... (6) 

in d > 3 dimensions [18] but one may also take into consideration the centrifugal-like spike in the 
potential (4). Thus, a generalization (1) — > (4) is to be understood as a transition to the singular 
models with F ^ and/or with 

L(L+1) + G = £(£+1) ^0, £=J G +(l+1) -J. (7) 



d L(L + 1) _ _ , , 
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In ref. [14], a theoretical merit of such a step has been seen in the identification of the older regular 
model (1) with the mere special case of eq. (4). Indeed, the vanishing of the charge F = F^ ES "> as 
postulated in ref. [1] results, in fact, directly from the QES conditions at £(£ + 1) = 0. 

Another consequence of the formal presence of the centrifugal term in eq. (5) lies in the related 
possibility of a modification of the potential (i.e., of the dynamics) by a mere formal change of the 
variables in eq. (5) [20]. This idea will not be discussed here in any detail but the interested reader 
may consult ref. [22] for an illustration. 

2.2 |x| 3> 1 asymptotics for the decreasing quartic potentials 

The general QES recipe starting from a polynomial potential [say, (1) or (4)] constructs its QES 
bound states [exemplified here by eq. (2)] in a way described by Magyari [23]. Basically, the 
construction parallels the harmonic-oscillator factorization tp(x) ~ cxp(— x 2 /2) x a polynomial 
where, for the VT— symmetric quartic model (4) with five real couplings, one extracts and separates 
the \x\ 3> 1 asymptotically dominant part of the normalizable (i.e., bound-state) wave function into 
its exponential factor, 

^(x)=cxp(-Ux 3 -±[3x 2 -nx) f>„(iz)"+P (3 = B/2, 7 = (J3 2 - C)/2 . (8) 

^ ' n=0 

In such a scenario and in a way extending eq. (2) to £ ^ 0, all the QES states will be characterized 
by the exact reduction of the infinite series to a polynomial, 

UJN+l = u N+2 = . . . = . (9) 

This means that our Schrodinger eq. (5) must be integrated over a complex contour of coordinates 
ieC which is bent downwards, say, towards its asymptotes 

Cleft ~ -ge+'W, Cnght^+ge-^', 0<<PUft,riglU<^ (10) 

with the large and positive real parameter g — > oo. Indeed, it is easy to verify that the exponent in 
eq. (8) decreases along both these half-lines, 

e'i ix ift = e +iie 3 (co S 3 W +isin3 W4 ) ^ e -i e 3 sin W x fl bounded oscMa tory factor, (11) 

e - I ' x right = e ~ 3 i<? 3 (cos3 Vrlght -i sin3 VrI9M ) _ g - i e 3 sin v „ ght x fl bounded osc ill ato ry factor. (12) 

As long as our problem is analytic in the whole cut complex plane (with the cut starting in the 
origin and oriented upwards), the contour C may be chosen as safely avoiding the singularity in 
the origin. Thus, in terms of the effective angular momentum £ we have to fix the sub-exponential 
exponent p = — £ in eq. (8) in a way compatible with both refs. [1] and [14]. 

2.3 Magyari's QES conditions 

The insertion of the ansatz (8) + (9) in Schrodinger equation (5) fixes the QES-compatiblc value 
of the coupling at the linear potential term, 

D = D(N) = 2{£ + (3 1 -N -I) (13) 
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and imposes, furthermore, the following overcomplete linear algebraic set of N + 2 constraints upon 
the N + 1 (arbitrarily normalized) wave function coefficients ui n , 

(2e-n){n+l)Lu n+1 + [F-2j(£-n)]Lu n + [E-j 2 +f3(2£-2n + l)}Lu n - 1 + 2(N + 2-n)Lu n - 2 = (14) 

where n runs from till N — 1. With obvious abbreviations, these equations may be re- written as 
a non-square matrix problem 



( S (F) U 

Ti(E) Si(F) Ui 

W 2 T 2 {E) S 2 (F) U 2 

W 3 T 3 (E) S 3 (F) 



U 3 



W N T N {E) S N {F) 
V W N+1 T N+1 (E) J 

It must be solved numerically in general [24]. 



UJ 
U>1 



= 0. 



(15) 



3 Two feasible versions of the QES construction 

3.1 The domain of the small G, L and £ 

In the practical computations one may treat eq. (15) as the two linear square-matrix eigenvalue 
problems 

( 5*0(0) U \ 

T^E) 5i(0) Ui / u \ / w \ 

W 2 T 2 (E) S 2 (0) U 2 ,. ( a.-j 



V 



T n ^(E) SW_i(0) U n ^ 
W N T N (E) S N (0) J 



= -Fe{E) 



\ UJ N J 



(16) 



( T!(0) S 1 (F) U x 
W 2 T 2 (0) S 2 {F) 



U 2 



V 



Tjv_i(0) S N - X {F) U N - X 
W N T N (0) S N (F) 

W N+1 T N+1 (0) ) 



( UJQ \ 

U>1 



\ UN J 



-E e (F) 



( LJQ \ 

U>i 



which are complemented by the two mutual-coupling conditions 

E = E e (F), F — F e (E) . 



(17) 



(18) 



An important note is to be added, based on the inspection of eq. (14). It reveals that in the 
Magyari's non-square matrix constraint (15), the upper-diagonal coefficient U n = {21 — n)(n + 1) 
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can in fact vanish at n = n(£) = 2£, i.e., at all the half-integer effective angular momenta £. In eq. (5) 
this represents a constraint upon the freedom in the choice of the spike strength G. Thus, whenever 
the integer 2£ does not exceed the dimension N, at least one of the coupled secular determinants 
factorizes [14] and a significant reduction of the complexity of the algebraic QES conditions is 
achieved. Still, the difficulties with the construction grow fairly quickly with the growth of 21 even 
at the integer values of this parameter. 

In the light of the latter comment one may be quite surprised by our forthcoming present main 
result saying that a dramatic and drastic simplification of the recipe recurs in the asymptotic domain 
where £ — > oo. 



3.2 Quartic QES models and their unexpected duality at the large I 



A few non-numerical samples of the solution of eqs. (15) may be found in refs. [1] (using £ = 
F = G = and several small N) and [14] (using £ = 1/2 and N up to 4, or I = 1 and N up to 
3). Also the results of these studies confirm that serious computational difficulties arise and grow 
very quickly whenever £ grows beyond one. In parallel [24], the form of cq. (15) appears to be 
perceivably simpler whenever all the values of the other parameters N, (3 and 7 become negligible 
in comparison with the partial- wave index L and/or with the strength of the core G. In the latter 
dynamical regime where I » max(7V, |/?|, I7Q) we may omit the negligible terms from our eq. (15) 
and get the leading-order version of the QES requirement, 

/ F-2 7 £ If. \ 



E + 2(3£ 
2N 



21 
F-2j£ 



4£ 



E + 2(il F — 2~/£ 



UJ 
U>1 



0. 



(19) 



\ WAT / 



6 E + 2f3£ F-2j£ 2N£ 
4 E + 2f3£ F — 2~/£ 
\ 2 E + 2[3£ J 

In this equation we may re-scale the coefficients u> n = h n £~ n / 3 and subtract the leading-order 
asymptotic approximants, 

F = 2j£ + 2s£ 2/3 , E = -2f3£+2t£ 1/3 . 



(20) 



This replaces eq. (19) by its strictly equivalent but strikingly simpler form 



/ s 
t 

N 



\ 



V 



N 
s 
t 



( h a \ 

hi 



0. 



(21) 



The phenomenologically most important and formally most remarkable consequence of this result 
is that it represents another manifestation of the Buslaev's and Grecchi's [18] duality between 
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Hcrmitian and non-Hcrmitian quartic oscillators, this time on the level of their respective QES 
subsets. Indeed, in the light of ref. [25], the same equation played the same role for the Hermitian 
asymptotically growing potentials 

V {Hermltlan) (r) = +r 4 + Br 3 + Cr 2 +Dr + F/r + G/r 2 , r G (0, oo) . (22) 

Amazingly enough, all the numerous differences between the potentials (4) and (22) (the latter 
being defined on the half- axis of course) disappear on the level of constraint (21). This enables to 
make the rest of our present text short. We may just cite the final (though, by the way, not so 
easily derived!) results of the extensive computations as performed in ref. [25]. In particular, this 
enables us to summarize that the real roots s = t = t(N) of cq. (21) form the N— dependent and 
equidistant multiplets of integers, 



t(N) =t k (N) =N-3k, fc = 0,l,.. 



(23) 



This means that the physically acceptable solutions of our present VT— symmetric t ^> 1 QES 
problem exist and occur in the multiplets with the following asymptotic energies and charges, 



E = -2(3£ + 2(N-3k)£ 1 / 3 + ..., F = 2j£ + 2(N - 3fc)^ 2/3 + . . . , k = 0, 1, . . . , 



• (24) 



One may now return to the elementary recurrences (21) and evaluate, very easily, the coefficients ui n 
of the wave functions in the same next-to-leading order approximation. In the light of the existing 
thorough analysis of this problem in the dual Hermitian context [25] , this task may already be left 
to the readers as an exercise. 



4 Summary 

Could we view the quasi-exactly solvable VT— symmetric quartic potentials as a choice, in some 
sense, "next" to the popular harmonic oscillator? In our paper wc tried to support an affirmative 
answer. 

During our study we felt particularly motivated by the technical difficulties arising in connection 
with the explicit construction of the quartic QES charges. Although, implicitly, they are defined by 
the coupled pair of the Magyari's polynomial algebraic equations for two unknowns, their practical 
determination must usually rely upon the computerized, Grobner-basis-based algebraic manipu- 
lation techniques and numerical root-searching [26]. In addition, it is quite unpleasant that the 
complexity of the latter algorithm grows fairly quickly with the growth of the degree N of the 
polynomial wave functions as well as with the growth of the angular momentum I. Finally, the 
situation significantly worsens whenever i ceases to be a half- integer [14]. 

We were encouraged by the well known fact that, quite often, the dependence of bound state on 
the angular momentum may get simplified in an asymptotic regime [27]. In the latter direction, our 
attempt proved successful. We found that several large— £ properties of our non-Hermitian model 
(like, e.g., the subtle QES-related cancellations of the separate elements in the infinite power series 
in x) find in fact quite close parallels in its self-adjoint predecessors. Many differences (e.g., the 
occurrence of complex couplings or the deformability of the integration contours in VT— symmetric 
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case) proved inessential. We arrived at the final version (21) of the Magyari's equation which is, 
from the formal point of view, identical with the equations encountered in Hermitian cases (so, we 
could also employ their well known solutions in our construction). 

In conclusion we may add that whenever necessary, one may leave the asymptotic domain and 
switch attention to the finite effective angular momenta I -C oo. The necessary mathematics may be 
found in the modified Raylcigh-Schrodinger perturbation recipe as adapted to non-square matrices 
in ref. [28] . It is worth emphasizing that it makes the full use of the finite-dimensional character of 
the Magyari's re-formulation of our present negative-quartic QES Schrodinger equation. For this 
reason it may be recommended as an efficient and systematic source of higher-order corrections. 
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